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ABSTRACT 

A new type of frame saw machine with four-bar parallelogram linkage mechanism has been recently studied 
around the world. In this paper, an advanced model including 8 control parameters, 9 functional constraints and 9 objec- 
tive junctions was developed with a goal to interact and deal with inconsistencies among four stages within design 
process of this machine, such as concept, analysis, technology and multi-criteria decision-making processes. The last 
stage or the most important one was in fact a multi-objective optimiuition problem; in order to tackle with it the authors 
used a visual interactive analysis method (VIAM) with an application of single-objective optimhation techniques. By 
using VIAM for tliree manujacturing scenarios, 28 Pareto optimal solutions have been determined; hence the designer is 
able to make a proper decision at every scenario. 
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INTRODUCTION 

A new type of frame saw machine with four-bar parallelogram linkage mechanism has been recently 
studied and developed [1-3]; an illustration of which is shown in Figure 1. Unlike traditional frame saw machine, it 
was designed to ensure dynamic self-balancing system; consequently, it can work steadily at high speed 
(approximately 3000 rpm) without consideration of an extra balancing mechanism. The movement is transmitted 
from the leading shaft (lower) to the shaft being led (upper) directly by mean of the saw blade without using any 
other drive system. The operational principle and structural improvements have significantly reduced the sawblade 
length in comparison with that of the equivalent traditional frame saw machine, thus stiffness and stability of the 
sawblade have been enhanced considerably. 

However, a task to solve inconsistency issues arisen within a new type of frame saw machine design 
process seems to be a complex multi-objective model. If the current multi-objective optimization methods are 
applied to deal with this, it might make the task becoming out of its essence. The reason is that, most of these 
methods implement the idea of converting the criteria into an equivalent function (or scalar methods) by using 
many techniques such as weighted minimax (maximin), compromise programming, weighted sum, bounded 
objective function, modified Tchebycheff, weighted product, exponential weighted sum, etc. [4-8]. Generally 
considering, the equivalent function of all criteria or objective functions is as 
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.. *Wa) = coMV 0 /Mft'o>!{<I> , a = ... <1> 

iollows: v ’ 1 v ' ’ v n , where v - a vector oi N parameters, ' - an objective 

function i (i = 1.. M). Afterward, it needs to find the extreme of the equivalent function; in order to do that, several current 

single-objective optimization algorithms can frequently be used [9-14], such as numerical direct methods or modern 

methods with behavioral characteristic-based approach to biological and molecular systems, swarm of insects, neural 

network etc. These algorithms yield optimization values with insignificant error. Nevertheless, by using the minimum of 

the resultant equivalent function, the outcome of any objective function from one algorithm is different from the one of the 

same function from another algorithm. In other words for example, the minimum of the resultant equivalent function 

\J/® ~ \J/® 5- ip® 

yielded from the algorithms are more or less equal, as assigned 1 ~ 1 ~ ' " ,e ' Wir respectively; while applied 

conversely, the results of every single-objective function from every algorithm are rarely alike, or 
(<J) :^(<|>.) @ ^ (<J>.) @ 

' 1 ' methodl ' '' method 2 " ’ ' 1 ' methodK | 



Figure 1: A New Type of Frame Saw Machine with 6 Saw-Modules 

Nonetheless, it is concerned that whether the results of every single-objective function are in coincidence with the 
requirements within the design process, its importance grade might be varied at a specific manufacturing circumstance. 
Besides, non-scalar optimization methods have been lately used, Brito and his co-workers [15] considered a relaxed 
projection method, which coordinates a retlection technique to obtain a feasible point with projected sub-gradient method; 
while in other study [16] the multi-objective optimization methodology based on harmony search and Pareto front 
approaches was proposed to design energy-efficient shading devices; and there are also other methods such as 
lexicographic, goal programming, physical programming, genetic, nash arbitration, described elsewhere in Ref. [17-18], 
and/or multi-objective genetic algorithm [0], Gaussian process regression [20], Nash/Adjoint optimization methods [22]. 

Nevertheless, most of the aforementioned methods did not concentrate on arranging an interactive panel, which 
allows the designer to be aware of feasible criteria. This leads to the fact that there is no basis to evaluate and analyze the 
resultant optimal solutions. Yet, these methods are only capable of finding out solutions in one-way direction; it means that 
although the resultant solution from every algorithm is Pareto solution, it might not satisfy the designer. For example, there 
is a case that one objective function yields a very good result, while another function has not yielded a desirable result yet. 
This shows that there is no professional interaction or control from the designer during solution search, that the obtained 
solution is barely the result from a fixed algorithm and no more beyond that. Moreover, these methods often tend to 
improve the efficiency of Pareto solution search (less search time, less amount of test points). It is noted that the most 
important thing is that the obtained solution needs to meet the technical requirements; and when they are varied, it must be 
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flexible to obtain a satisfactory solution correspondingly. Therefore, in order to deal with the complex multi-objective 
model within design process for a new type of frame saw machine, the authors propose a visual interactive analysis method 
or VIAM based on a single-objective optimization techniques. 

PROBLEM STATEMENT 

The main structure of the frame saw machine consists of similar sawblade modules, as shown in Figure 2. 
The sawblade module is installed on lower- and upper shafts, and moved by four-bar parallelogram linkage mechanism. 
There is a key connection between the shaft and the eccentric disc. Besides, the eccentric disc is joined with the housing 
part by hinge connection. The rotational movement of the upper shaft is yielded from the lower shaft with rotation speed n 
(revolution per minute) through the sawblade modules themselves. The matter-spot on saw-module moves with 
eccentricity of circular motion, e, at the same velocity and acceleration [23]. Centrifugal acceleration of the matter-spot on 
saw-module has a constant magnitude such as a=e*(2jt/60) 2 , but its direction changes continuously. During no-load 
operation, the sawblade is subjected to initial tensile force F 0 with eccentricity e 3 , inertial uniformly distributed force q 
along blade length, and moment generated by the inertial force F b of counterbalance. Lower - and upper housing parts are 
designed to ensure that their gravity centers coincide with the center of eccentric disc (A h A 2 ), thus they do not cause 
inertial moment acting on the sawblade. The entire moment of inertia effects on the blade with a maximum magnitude at 
the location, where the centrifugal acceleration is in the horizontal direction, i. e. the position with rotational angle a = 0° 
or 180°. Assumed that the machine needs to be designed to be able to saw logs with the largest radius h, the minimum free 
length of blade and the total blade length would be L = li + 2e and L 0 = L + 2L k respectively, where L k is the partial blade 
length, which is clamped at each end and is considered as a constant in the calculation model. 

The set of control parameters, which decides criteria of saw machine, includes: eccentricity of the circular motion 
e and shaft rotation speed n, blade dimension, counterbalance mass m b and distance h b , the magnitude of initial tensile 
force F 0 and eccentricity of saw blade tension e/, as described in Table 1. Besides, the functional constraints and quality 
criteria or objective functions are presented in Table 2 and Table 3 respectively [24]. 


Table 1: Control Parameters 


Parameters 
-Symbols in 
Model 

Original 

Symbol 

Initial Minimum 
Permissible Value 

Initial Maximum 
Permissible Value 

Units 

Definition 

0-1 

E 

0.03 

0.035 

m 

Eccentricity of the 
circular motion 

0-2 

b 

0.06 

0.1 

m 

Saw blade width 

a 3 

t 

0.001 

0.002 

m 

Saw blade thick- 

ness 

a 4 

ei 

0 

0.08 

m 

Eccentricity of 
saw blade tension 

a 5 

h b 

0.1 

0.2 

m 

Distance h b 

a 6 

m b 

0 

1 

kg 

Counterbalance 

mass 

a 7 

F 0 

500 

2000 

N 

Tension force 
magnitude 

a 8 

n 

2000 

3000 

rev./min. 

Shaft rotation 
speed 


www.tiprc.ors 


SCOPUS Indexed Journal 


editor@tjprc. org 

















452 


Dang Hoang Minh, Phung Van Binh & Nguyen Viet Duc 




a) Stmcture sketch 

1,9- Counterbalance; 2 - Upper Eccentric; 3 - Upper Shaft; 4 - Upper Housing part; 5 - Saw Blade; 
6 - Lower Housing Part; 7 - lower Shaft; 8 - Lower Eccentric 


Figure 2: Saw Blade Module in a New Type of Frame Saw Machine 


Table 2: Functional Constraints 


Symbol 

Requirement 

Explanation 

fi 

<0 

Requirement for an absence of resonance 

f 2 

>0 

Stability requirement of sawblade under inertial forces 

h 

<0 

Balance state of the saw-module 

U 

>0 

Tensile force requirement 

fs 

>0 

Strength requirement of sawblade 

f 6 

>0 

Fatigue requirement of sawblade 

h 

>0 

Initial rigidity requirement of sawblade 

f* 

>0 

Stability requirement of sawing processes 

f 9 

>0 

Specified requirement (real value exists) for criteria <f>3 


The mathematical model herein includes 8 control parameters (Table 1), 9 functional constraints (Table 2) and 9 
objective functions (Table 3). It is necessary to find solutions satisfying with the constraints and the functions need to be 
optimized in compliance with the technical requirements of the saw machine. 


Table 3: Objective Functions 


Symbol 

The Direction of 
Improvement 

Explanation 

cDj 

MIN 

Total mass of the sawblade and the 
counterbalances 

(J >2 

MIN 

Overall dimension 

O 3 

MAX 

First natural frequency of the saw-module 

O 4 

MAX 

Critical speed of shaft rotation 

®5 

MIN 

Tension force magnitude 

®6 

MAX 

Operating speed of shaft rotation 


MAX 

Stability of sawing processes 

tpg 

MAX 

Initial rigidity of unstrained sawblade (in the 
absence of sawing force) 

O 9 

MIN 

Sawblade thickness 
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_ _l 


(a) Trend I (B) Trend II 

Figure 3: Algorithm of VIAM 

METHOD AND ALGORITHM OF SOLUTION 

The main idea of VIAM is to use the single-objective optimization techniques as a tool, with the aim of finding 
suitable solution for multi-objective model, satisfying designer requirements. The detailed steps in solution search and 
generalized diagram are described in Figure 3a and 3b, respectively. 

Function Value and Interactive Panel 

The model { 1} is the first step, as shown in Figure 3a, hence it is necessary to optimize the vector <I> of M 

objective functions with the parameter vector ^ , as well as the functional constraints ^ (where, a = {a h a 2 ,..., a N ), 
f = f(a) = { /i(a),/ 2 (a),...,/x(a)}, <D = <t> (a) = { <T>i(a), <I> 2 (a),..., ® M (a)}). The condition of parameter a is as follows: 

a ; < a t < li . Oj - a jQX a t < b k ^ ^ £ [l; }V]). Besides, the condition of functional constraints f is as follows:/(a) < 0;/„(a) 
= 0;/,(a) < 0 or/ 0 (a) / 0 (/, m, n, o e [1; A])- Next step {2} is to use up-to-date single-objective optimization algorithms 
[4-9, 12, 13], with the aim to determine minimum and maximum values or MIN<t>, and MAX<f>„ respectively, of an 
objective function (taking into account the relationship), and include these values into a table. This is an interactive panel, 
which helps engineers significantly to analyses and makes a decision during solution search. When the final solution is a 
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<D s ={<frf, *?,...,*£} , . < f ) f’ e [MIN<5.;MAX<5.] 

vector 1 ‘ with values satisiying designer requirements, 1 . e. , there are two 

principal trends in solution search of multi-objective optimization problem. 

• Trend I (Figure 3a): Based on the table in the step {2}, the designer defines a “strict” priority procedure for 
criteria, for instance in {3}: (the rest criteria (M - P) are not considered, because their importance grades do not 
reach to a vigilant level). It means that the function 1 is the most important one, it needs to be the most optimal 
value, then the function 2 also has to reach the most optimal value after making concessions on importance to the 
function 1 , and the process carries on for other subsequent criteria. 

• Trend II (Figure 3b) is equalizing classification, i. e. there is a group of first priority functions, which have a 
similar importance, such as {12}: Dealing with an optimization of one function in this group might aggravate the 
results of other functions, thus in this circumstance the analysis and discussion among designers need to be carried 
out. 


It is important to indicate that the Trend II is a generalized case; it can be converted into the Trend I, when Q = 1 
occurs gradually M times, or extend to all of functions as Q = M. However, in many cases though Trend II is used to solve 
the problem, but when it is necessary to make a decision, sometimes randomly it is prone to Trend I, because it is rarely 
capable of optimizing two or more functions simultaneously, consequently there needs to be a “reluctant’’ agreement to 
make a priority for any individual function. 

Block “Approach by Priority Order of the Criteria” 


This approach to finding the solution by Trend I (Figure 3a) is as follows. First, it needs to optimize the most 
important function 1 (step {4}). Hence, the optimal value of function 1 or MIN<T> can be used for further steps. However, 
while setting a goal for <t>, is such high, the obtained solution of other functions might be undesirable. Therefore, the 


1<I> 1 = MIN<I> +AE 

designer needs to extend the requirement for <t>], and set its value within a threshold L 1J 1 1 (with the 

AE >0 <T>. < [<!>, 1 

deviation 1 ). If so, among the obtained solutions that L J , there might be the most favorable option for the 

rmri$ $ — [(J) ]<n 

function <t> 2 . While searching the optimal value for the function <t > 2 or 2 , it needs to add the constraint 1 L 1J 


(an absolute form) or 


[(*i - [*i])/|0-5-■(«!*. +[<h,])|] <0 


(a relative form) {5}. As soon as at this step, the function 2 is opti- 


mized with consideration of the constraints of the function 1 , thus the obtained value nlm< ^ 2 is often greater than the 

MINcE> min^ 

2 in the interactive panel {2}. Then, the value 2 needs to be verified by the experts, because there might be a 

[q> 1 

situation that this value is inferior and inacceptable. If so, the threshold L IJ needs to be extended more with the aim to 
re-solve the step {5}. While, in the case that the obtained value min<I >2 complies with the requirements, similarly they set a 

threshold m ' nC b 2 + ^ ^ f un ction 2 (with the deviation y ®) with the goal to find solutions for optimizing 


the succeeding functions (step { 6 }). While searching the optimal value for the function 3 or 


min<I>, 


, it needs to add two 


constraints 


|(<f>, -[«>,] < 0) |[(4>, -[«>,])/|0.5 ■ (4>, + [4>,])|] < 0 

[(<h 2 -[<f> 2 ]< 0 ) 01 |[(<f > 2 -[<f> 2 ])/| 0 . 5 -(<f > 2 + [<f > 2 ])|]<0 (step{7}) 
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r<j> ] = min<I> + \F. 

By the same manner, the designer sets a threshold L 3J 3 3 for this function (step {8}) etc., and so on 


there are a series of thresholds 




for the first (P- 1) functions (step {9}). Based on ( P-\ ) constraints of the 


corresponding threshold and necessity of minimizing P" function, eventually the solution with a desirable priority order 

4>® = {<J>f, 4>f,...,<l>®} 

can be achieved (step {10}). The obtained result L ' is just a Pareto solution ({11}), which means that 

there is no better solution for all functions simultaneously. In case, the designer needs to verify exactly whether the 
obtained solution is Pareto one, it is possible to shift to Trend II at step {23}, which will be described with more details in 
section 3.3. 


Note: When dealing with the problem by this approach, the most important issue is a next question. Although an 
individual function might reach an optimal value, whether or not this value is desirable, while extending slightly the 
threshold, it helps to acquire more optimal solutions for other criteria? This cannot be predicted by an engineer; therefore it 

fd>.l = mi n + \K 

is necessary to have a tool to control the threshold L ' J ' ', so that it is capable of extending constraint condition 

of subsequent functions relevantly and precisely. Thus, it is able to obtain suitable solutions for other functions. 

Block “Approach by Equalizing Classification of the Criteria” 

Trend II (Figure 3b) suggests that there is a group of Q functions with the same importance, whereas the rest func- 
tions are random (step {12}). The experts herein need to set a “threshold” value of [<t>,] for the function i in the interactive 

OX = {<I>X 1 ,<I>X 2 ,...,<I>X 0 j OX 

panel {13}. Range of satisfied solutions will be a set of vectors ' , where the final value 'of the 

<hX e min<I> •[<!> 1] 

function needs to comply with the condition ' L J . The main idea is to find that vector ©X. This approach 

presents as follows; firstly considering the aforementioned condition as the constraints {14}, hence it is necessary to find 
not only N parameters Xj (j = 1.. N ), but also Q final value of the function ©X, (i = 1.. Q) while minimizing the equivalent 
error function 


min F = 


F: 


-<hX.)/ 0 . 5 (q>, +4>X,)| 


min £[(<*>, -® x /)/ 0 - 5 (®i +4> x /)] 


£(+-<I>X, 


relative form {14} or 

A ={ a ,©X 


min F = min -j 'j' |<f>, - <1>X,. 


absolute form. In other words, the vector parameter is 


(W+(2)xl [Aixl <2x1 


together with its 




condition. The aim of using the final value of the function ©X, instead of the threshold L ' J in the equivalent error function 
F is to find the solutions being better than those set by the threshold (if they exist). 


In case, the optimal value of the function F approaches to zero (step {15}), this means that there might be a 
presence of several satisfied solutions. Here, either the threshold value [©,] is possibly still broad or there is a random 
appearance of one solution. There are two ways to deal with this circumstance. The first one is to consider the threshold as 
Q constraints, and then acquire Q Pareto optimal solutions by determining the extreme of each function (steps {16} and 
{17}). The second one is to use condition of function F 0 minimization, which is an equivalent penalty function, in order to 
find the satisfied solution at the steps {18} and {19} (this technique is described in details in the Section 3.5). Among the 
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satisfied solutions, it needs to apply a “filter” to find the most optimal one by Pareto. Taking into account that the filtered 
solutions might not be the global optimal Pareto solutions, but they relieve the designer and they are the best existing 
solutions. If there is an available time, it is possible to verify whether they are the global optimal Pareto solutions (step 
{23}) by considering them as a threshold and return to the steps {13}, {14}, {15}. 

Nevertheless, while dealing with this type of problem in practice, the aforementioned favorable circumstance 
rarely appears, and the typical circumstance is when the minimum value of the function F could not approach to zero. This 

is a case of step {20}, min F = E°>e, j t corresponds to the parameter vector a and the final values of the function . 
Once it is perceived that the up-to-date single-objective optimization technique is strong enough to determine min/ j and 
min/- does not approach to zero, it means that within the threshold value, indicated by the designer, there is not any 
satisfied solution. Here, there are three following questions: 


• First, whether the threshold value of the function [<t>i], which is included in the table {13}, is appropriate or not? 

• If truly there no exists any satisfied solution, what to do next? 

• If it needs to alter the threshold value of the function, of which one needs to do? How does it alter (increase or 
decrease)? And at which magnitude? 

Generally speaking, the designer set the threshold value mostly based on their practical expertise, thus upon 
defining the solution, they themselves do not know the threshold indicated by them whether achievable. Thus, the step 
{20} helps them to know the thresholds indicated by them being unachievable simultaneously. After checking the devia- 
e =<f>.(a°)-<f>X° <f> (a°) 

tion 1 v ' ' among the values ol the function ' v m comparison with the final value of the corresponding 

<f>X° <f> (n°) 

function ' (step {21}), the designer will have two options; First, the functions ' v do not coincide with the prelimi- 

(a =a°,<f> (a°)) 

nary observation, but they satisfy the designer, and consequently the obtained solutions v ' are acceptable. 

Or second, the designer needs to alter the threshold of the function in order to find more suitable solutions. In case it needs 

£ 

to alter the threshold, it will be done by using the deviation ' obtained at the step {21}. The negative or positive value of 
£ 

the deviation ' will decide how to alter the threshold; if it is positive, it is required to increase, and inversely to decrease. 

I £ 

The value 1 : decides the magnitude of alteration. If the deviation ' is approximately zero, the threshold of the function 
in relation to that of other functions is suitable and there is no need to alter. The preliminary threshold turns to be 
[<f> 1 =[<f> ) + £■ 

L L ' J '({22}), and it needs to be included into the table {13} in order to repeat a calculation process again. 
Accordingly, when the new threshold satisfies the condition {15}, it is necessary to carry on the next steps {16} and {18} 
to find the satisfied and Pareto solutions. It is important to note that the alteration of threshold at step {22} is crucial; it 
decides whether the satisfied solution exists. 


Spatial Parameter Conversion Technique 


6 = £--£5 

In order to avoid the difficulty of optimization algorithms, when the parameters ( ’ vv< -' ,xl " Vxl 0x1 * , (i = 1.. Q) are 
set in a limited value domain, YIAM uses a spatial parameter conversion technique at the steps {1}, {5}, {7}, {10}, {14}, 
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{16}, described in details elsewhere [14]. In Figure 4, the continuous line presents the acceptable value domain of 
parameter, while the dashed line illustrates the unsatisfied value domain. The spatial parameter a (limited) is converted 
into the spatial parameter t (unlimited) in three cases, each of them may use a different conversion method. 


(a) 


oc 


t 




i b k / 
s 

t • — 




a. 



Figure 4: Spatial Parameter Conversion Technique. 

(a) - Limited Spatial Parameter a; (b)- Unlimited Spatial Parameter t 


{+ 

,0x1 

T =' 


[lVxl 


is referred as ( ' v7e,xl 

[Axl Qxl J 


®T = {<UT„<1>T 2 ,...,<1>T 0 } 

/ 2 ,..., /,v}, - . In the new space T, the parameter is varied from -co to +co, hence it is not limited as 

before in the space A. This helps to improve the stability in all optimization algorithms used, because there might be any 
invalid solution that makes the search process stop or error occurs. 

Technique by using a Penalty-Function to Determine Pareto Solutions 

The penalty-function technique is a proper approach to fmd the solutions at the steps {18} and {19} in VIAM. 
Range of satisfied solutions is determined by using a condition of function F 0 (an equivalent penalty function) minimiza- 


i\ a 

F 0 (A) = y) Penalty ( a ..) + ^ Penalty (/, (n)) + V} Penally 


tion. This function F 0 has a generalized type as 


0-5>,(a) + [*,]| 


with a 


spatial parameter A, or 


F 0 (T) = y Penalty (/, (t)) + J^Penalty 

i =1 1=1 


0-5'| < h / (t) + [ < f> / ]| 


with a spatial parameter T. In order to find 


many satisfied solutions at the step {19} with an application of single-objective optimization techniques, it is necessary to 
establish a sufficient small increment in search process and to use a significant number of initial test vectors [5]. 


It is noteworthy that the options based on the use of an equivalent error function and of an equivalent 
penalty-function F 0 at the step {18} are nearly similar, as they both yield the satisfied solutions. Therefore, it is possible to 
use flexibly one or another option based on a particular circumstance. However, the condition of an equivalent error 

function mm F - £ is not as strict as the one mm _ **. Consequently, at the step {14} if m ' n ^ - £ , this is a signal, which 
indicates that there is a high possibility of satisfied solution existence, and the step {18} needs to be used many times with 
numerous initial test in order to determine potential solutions. 


RESULTS AND DISCUSSIONS 

Preliminary Step 

Determine the maximum and minimum values of each function and include them into an interactive panel {3}. 

By using a single-objective optimization method, the extreme of 9 functions can be found, and they are included 
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into the interactive panel Table 4 below. It is noted that the functions 3, 4, 6, 7, 8 need to be maximized. However, to 
facilitate the analysis, it is prone to find the minimum value of argument. Also, range of 9 functions has been determined 
and provided in the corresponding column of the table. Based on that, the designer set a required threshold for every 
function. In principal, the function would approach to the minimum value, which is presented in the upper row of the 
Table 4. The minimum value of the function 4 is theoretically -oo, though here with the actual calculation tool it is consi- 
dered as -2.89E + 78. This means that critical speed of shaft rotation may proceed to extremely enormous. 

Scenario 1 - Trend I 

Based on the analysis and discussion, the designer make a decision on this particular scenario about priority order 
of the criteria as follows: 

{4>, h+ 4> 6 h+ 4> 5 h+ <b 9 h+ 4>, h+ 4> 4 h+ 4> 7 h+ 4> s h+ 4>,} 

Solution search procedure according to Scenario 1 is included in Table 5. The highlighted green color cell 
indicates the “good/best” obtained value of the criterion, which is being optimized deliberately. While, the highlighted red 
color cell implies the “bad/worst” obtained values of the rest criteria (if any). Eventually, the highlighted purple color cell 
represents the “relatively good / valid” obtained values, which are accepted by the experts (if any). 


Table 4: Interactive Panel {3} 


MIN©! 

min © 2 

min © 3 

MIN©4 

MIN©5 

min © 6 

min © 7 

min © 8 

min © 9 

0.371 

0.855 

-124.85 

-2.89E+78 

871.2 

-3000 

-4003.8 

-122700 

0.0011 










[®il 

[©2] 

[© 3 ] 

r©4i 

r© 5 i 

[©6l 

[© 7 ] 

[©sl 

r©oi 










MAX©i 

MAX©2 

MAX©3 

MAX©4 

MAX©5 

max © 6 

MAX©7 

max © 8 

max © 9 

2.7254 

1.065 

-50.12 

-3000.2 

2000 

-2000 

-997 

-49996 

0.002 

Mass 
-> min 

Sawblade 
overall 
dimension 
-> min 

Argument 
of “First 
natural fre- 
quency of 
the saw 
module” 

-> min 

Argument 
of “Critical 
speed of 
shaft rota- 
tion” 

-> min 

Tension 

force 

magnitude 
-> min 

Argument 
of “Operat- 
ing speed 
of shaft 
rotation” 

-> min 

Argument 
of “Stabil- 
ity of saw- 
ing 

processes” 
-> min 

Argument 
of “Initial 
rigidity of 
unstrained 
sawblade” 
-> min 

Sawblade 
thickness 
-> min 


In the Table 4, the optimal value of the function 1 is 0.371, but the designer observe that the value closest to 0.371 
can be accepted being 0.5. If so, it is possible to obtain more solutions, allowing for optimizing other functions. In order to 

[<j> 1 = 0.5 

verify whether there exists the valid solution complying with the threshold L IJ , it is necessary to use the 


minF = minjldi, - f<f> l l|} —> 0 

minimization conditions of the error function F or 11 IJ and one of the solutions (Sl) is turned out, 

as shown in Table 5. Though the solution (Sl) is not acceptable due to the fact that the functions 5 and 6 approach to the 


worst value, but the existence of several valid solutions at this step points out that the range of threshold value 
can be used for the next search process. 


[<&,] = 0-5 
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Table 5: Solution Search Procedure According to Scenario 1 


<D 1 = 

0.5 

1 

(Sl) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

o-i 

0.030642 

0.060499 

0.001635 

0.003476 

0.16386 

0.073975 

1963.6 

2116.8 

- 

fi 

-2541.2 

565.7 

-18.804 

1007.1 

29.337 

37.542 

8500 

723.3 

2842.8 

<D, 

0.5 

0.984 

-116.45 

-4023.7 

1963.6 

-2116.8 

-1223.3 

-58500 

0.001635 

®6 = 
-2900 

1 

(S2) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

a.j 

0.030656 

0.060012 

0.0013054 

0.01849 

0.2 

0.11071 

2000 

2963.7 

- 

fi 

-1358.7 

0.3 

-87.672 

1239.2 

5.6049 

1.3817 

104. 

509.9 

979.33 

cD, 

0.50024 

1.0563 

-108.06 

-4446 

2000 

-2963.7 

-1009.9 

-50104 

0.0013054 

<d 5 = 

1550 

1 

(S3) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

a i 

0.030001 

0.060001 

0.0016 

0.000008 

0.1999 

0.080006 

1533.9 

2900 

- 

fi 

-1481.6 

- 0.1 

-9.3543 

590.10 

20.153 

7.831 

-4 

547.5 

2235 

<D, 

0.50071 

1.0550 

-109.54 

-4349.9 

1533.9 

-2900 

-1047.5 

-49996 

0.0016 

®6 = 
-2800 

1 

(S4) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

°-i 

0.03001 

0.06 

0.0015971 

0 

0.2 

0.07993 

1550 

2900 

- 

fi 

-1491.2 

0.9 

-88.325 

692.96 

20.059 

7.806 

4 

545.7 

2228.8 

cD, 

0.49996 

1.055 

-109.78 

-4351.4 

1550 

-2900 

-1045.7 

-50004 

0.0015971 

®6 = 
-2760 

1 

(S5) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

a.j 

0.030009 

0.060005 

0.002 

0.00019 

0.19947 

0.036986 

1543 

2804 

- 

.f, 

-1952 

1.1 

-136.41 

553 

47.712 

58.548 

29362 

1290.6 

6430.6 

<D, 

0.49991 

1.054 

-118.9 

-4207.6 

1543 

-2804 

-1790.6 

-79362 

0.002 

O 9 = 
0.0017 

1 

(S 6 ) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

°-i 

0.03 

0.06 

0.001676 

0.00026 

0.2 

0.071454 

1550 

2774 

- 

f, 

-1673.2 

67.3 

-91.463 

666.92 

24.963 

12.248 

5002 

671.3 

2771.6 

<D, 

0.4998 

1.055 

-111.18 

-4262 

1550 

-2774 

-1171.3 

-55002 

0.001676 

O 3 = 

_ 111 <D 4 = 
-4200 

1 

(S7) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

OLj 

0.03 

0.06 

0.0017 

0.0001 

0.2 

0.068918 

1550 

2775.6 

- 

f, 

-1690 

56.5 

-94.426 

659 

26.242 

13.373 

6498 

708.8 

2958.7 

<D, 

0.49984 

1.055 

-111.65 

-4248.1 

1550 

-2775.6 

-1208.8 

-56498 

0.0017 

O 7 = 

-1237 

1 

(S 8 ) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

aj 

0.03 

0.06 

0.0017002 

0.00254 

0.2 

0.069008 

1550.1 

2788.4 

- 

fi 

-1676.4 

11.3 

-99.054 

659.03 

23.951 

12.799 

7515 

737.9 

2958.5 

cD, 

0.50006 

1.055 

-111.62 

-4199.5 

1550.1 

-2788.4 

-1237.9 

-57515 

0.0017002 

(S9) 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

a, 

0.03 

0.060028 

0.0017007 

0.00244 

0.19957 

0.068942 

1550.3 

2798.3 

- 

fi 

-1666.5 

0.5 

-99.786 

659.07 

23.989 

12.749 

7515 

738 

2965 

<D, 

0.5002 

1.0541 

-111.62 

-4198.2 

1550.3 

-2798.3 

-1238 

-57515 

0.0017007 


(<t> — r<j> i<o) 

According to Trend I (step {5} in Figure 3a), the constraint v 1 L u ' can be added to find the smallest value 
for the function 6. The solution (S2) is obtained with the best value for the function 6 being -2963.7. For the functions 1 
and 9, the affirmative values are also obtained; but this solution is useless because three functions 5-7-8 approach to the 
worst value. Thus, it is necessary to “slacken” or extend the condition of the function 6. While, the designer estimate that it 

[o ] = -2900 

will be a success, as soon as the operating speed of shaft rotation reaches to 2900 rpm, thus the threshold L 6J is 


(<E> -f<E> ]<o) 

chosen. With the constraint v 6 L 6J ’ and minimization condition of <f> 5 , the solution (S3) is turned out. Here, 
although for the functions 1 and 6 the valid values appear (there is min<f> 5 = 1534), for the function 8 there exist the worst 

[o 1 = 1550 

value being -49996, which is unacceptable. Similarly, it needs to loosen the condition and opts for L 5J in order to 
find the extreme of <f> 9 . However, the solution (S4) is not affirmative in comparison with the (S3) (<f>g= -50004 is very 
low), hence once more it is not accepted. At this point, it is notably necessary to loosen the conditions of the function 6, 


and after analysis the designer decide to continue the search process of valid solutions with L 6j ~ . And indeed, it is 

rewarding at (S5), when the function 8 has improved with the value of -79362, but if looking into the function 9, the value 
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reaches to the least of 0.002. At this time, the designer is forced to agree with the fact that the value decreases as 

anc | jjjg t0 p |j m j t 0 f function 8 is not allowed to be greater than -55000. Here, there exists the solution (S6), 
which complies with the condition of functions 1- 5 - 6 - 8. Although <f> 9 = 0.001676 differentiates very much from the 

solution (S2), it is the best one in relation to other functions. Therefore, 0.0017 () p tet | f ()| . con sidering as a basis to 

continue finding valid solution for the function 3. In the obtained solution (S7), the best value for the function 3 is -111.65, 


thus 


[^3] = -Hl ; 


is taken into account. Still, the function <t> 4 is not optimized better than the value of -4248.1 in the solution 
[O l=-4200 

(S7). Accordingly, from here L 4J is opted for finding the minimum of function (f> 7 . 

In the solution (S8), the function 7 reaches to the optimal value of -1237.9, and other functions are also valid. 
[<l> 1 = —1237 

Next, L 7J should be used to find the minimum of cD 8 and <f> 2 , and the solution (S9) is turned out. Since there is a 

slight discrepancy between (S8) and (S9), evidently these two solutions cannot be optimized any more with the priority 


order 


{<f>, h-> <f> 6 h-> <f> 5 h-> <f> 9 h-> <f> 3 h-> <f> 4 h-> <f> 7 h-> <f> g h-> <f>,} 


, which has been set by the designer initially. Hence, 


[<f>„l = -57515 [<f>,l = -1.054 , . , , . . , . , . , T ■ „ ■ 

L J and L J are those two optimal solutions tn accordance with scenano 1. It can be seen visually m 

Figure 5 that value domain of parameter u.\ + a 7 for the scenario 1 is “narrow” circumstance, it is almost a single value: u.\ 

= 0.03; a 2 = 0.06; 0.3= 0.0017; a^~ 0.00249; a 5 = 0.2; ae= 0.069; a 7 = 1550.2; only in case of parameter a s there is a range 

ofvalues [2788.4; 2798.3]. 

0.03,0.03 


=J> > 


0.06^ 0.060028 


^0017002^0.0017007 


=> a : 


0.00 244 ^ 0.00254 


E>°4 


0.19957,0.2 


( UI68942 ^0.069008 


=>"t 


1550.1 1550J 




2788.4 „ 2798.3 


=> a » 


Figure 5: Intuitive Chart of Solution Distribution in Accordance with Parameters - Scenario 1 
Scenario 2 - Trend II (Figure 3b) 

During another manufacturing circumstance, the designer prioritize the improvement of 4 functions at the same 
& < l ) & < 3 ) &O c l 

time such as 1 1 5 6 9J . These functions have a similar importance; for the remaining five functions a random 

value can be obtained. Based on experience, the designer set the admissible threshold values for the four functions in table 
{13} as follows: [(DJ = 0.8; [<t> 5 ] = 1000; [<t> 6 ] = -2800; [<t> 9 ] = 0.00147. Here, it needs to add to the constraint 4 extra- 
conditions MIN<f>i < <f>Xi < [<t>i], where i = {1, 5, 6, 9} and check the extreme condition for the equivalent error function 


as 


min F = min 1 / Z [(<!>,- < f>X,.)/0.5(<f>,+<f>X,.)] 2 l<£ = 10 


. While, there exists the solution complying with this condi- 


tion, this means that there might be more valid solutions. In order to find the satisfied solution, it needs to use the condition 
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F o ( T )=X Pe ” a M/( t )) + X Penalt y 

Olohjil min F — 0 _ /=i j= 1,5, 6 ,9 

0 together with the equivalent penalty function 
From this, after filtration there are 12 Pareto optimal solutions, which are included in Table 6 . These solutions are 

represented visually by the chart in Figure 6 a. The functions are converted into a relative form 

^ <t>' - MINT, 

(unitless) MAX<f> ; MIN<f> ; , . _ \ gj _ jq.. 21, that they have the same scale in the chart. It shows that there is a 

similarity on magnitude from the functions 1, 2, 3, 5, 6 , 7, 8 , 9. 


o.5>,(t)+r<f>,i 


Table 6: Pareto Optimal Solutions in Scenario 2 


NO 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

(S10) 

a, 

0.030008 

0.07868 

0.00147 

0.07537 

0.19822 

0.19475 

871.24 

3000 

- 

®/ 

0.8 

1.0515 

-75.763 

-4854.1 

871.24 

-3000 

-1432.2 

-50033 

0.00147 

(Sll) 

a, 

0.033209 

0.07445 

0.0014622 

0.07512 

0.19941 

0.19796 

968.57 

2814.5 

- 

®/ 

0.78772 

1.0602 

-70.578 

-5112.7 

968.57 

-2814.5 

-1426.7 

-50102 

0.0014622 

(S12) 

a, 

0.030152 

0.066552 

0.001469 

0.070721 

0.1993 

0.19976 

907.44 

2983.9 

- 

<t>/ 

0.74672 

1.0539 

-75.251 

-2.89E78 

907.44 

-2983.9 

-1315 

-50071 

0.001469 

(S13) 

a, 

0.03 

0.06258 

0.001468 

0.060869 

0.19569 

0.18433 

981.89 

2801.8 

- 

®/ 

0.6947 

1.0464 

-79.108 

-2.89E78 

981.89 

-2801.8 

-1257.2 

-50236 

0.001468 

(S14) 

a, 

0.03 

0.06 

0.001468 

0.059957 

0.19975 

0.15397 

989.69 

3000 

- 

®/ 

0.62054 

1.0545 

-79.15 

-5478.1 

989.69 

-3000 

-1231.3 

-50042 

0.001468 

(S15) 

a, 

0.032574 

0.065882 

0.0014563 

0.07112 

0.2 

0.18723 

994.6 

2800.8 

- 

®/ 

0.71882 

1.0602 

-70.038 

-2.89E78 

994.6 

-2800.8 

-1338.7 

-50098 

0.0014563 

(S16) 

a, 

0.030409 

0.064067 

0.0014598 

0.06573 

0.2 

0.1652 

968.21 

2999.9 

- 

®/ 

0.66291 

1.0558 

-75.668 

-5387.4 

968.21 

-2999.9 

-1284.7 

-50140 

0.0014598 

(S17) 

a, 

0.031466 

0.074149 

0.001469 

0.07999 

0.1901 

0.19659 

938.46 

2806.2 

- 

®/ 

0.78224 

1.0381 

-70.468 

-4443 

938.46 

-2806.2 

-1461.6 

-51793 

0.001469 

(S18) 

a, 

0.031049 

0.064156 

0.0014517 

0.07102 

0.17356 

0.21831 

976.8 

2823 

- 

®/ 

0.76868 

1.0042 

-70.697 

-2.89E78 

976.8 

-2823 

-1316 

-50376 

0.0014517 

(S19) 

a, 

0.034431 

0.075785 

0.0014698 

0.07342 

0.2 

0.19791 

997.78 

2801.5 

- 

®/ 

0.79884 

1.0639 

-70.939 

-4599.9 

997.78 

-2801.5 

-1447.7 

-50018 

0.0014698 

(S20) 

a, 

0.030021 

0.060413 

0.0014528 

0.06572 

0.12284 

0.24242 

985.9 

2870.5 

- 

®/ 

0.79635 

0.90072 

-74.039 

-4545.8 

985.9 

-2870.5 

-1257.9 

-50216 

0.0014528 

(S21) 

a, 

0.030069 

0.067572 

0.0014668 

0.08 

0.2 

0.16641 

859.01 

2840.7 

- 

®/ 

0.68469 

1.0551 

-71.323 

-4403.5 

859.01 

-2840.7 

-1348.5 

-50129 

0.0014668 


Regarding the function 4, only 4 solutions (S12), (S13), (S15), (S18) are superior to those from the other 8 
functions. Thus, it indicates that to meet the condition of scenario 2, four solutions (S12), (S13), (S15), (S18) are available. 
Similarly, Figure 6 b represents valid solution distribution in accordance with parameters in scenario 2. It can be observed 
that the parameter domain of valid solutions in scenario 2 is quite large, which helps to find more solutions but the 
magnitude of function is not much different from 12 abovementioned solutions. 

Scenario 3 - Trend II 

Assumed that the designer could not know priority order among functions, yet all of them are considered equally 
due to the own significance for the saw machine. This means that 9 functions have an equal importance grade. If so, the 
designer have to set threshold values for all functions, such as [<t>!] = 1.1; [<t> 2 ] = 1.0; [<t> 3 ] = -80; [<t> 4 ] = -5000; [<t> 5 ] = 
1300; [<t> 6 ] = -2800; [ct> 7 ] = -1400; [ct> 8 ] = -55000; [<U 9 ] = 0.0015. Similarly, it needs to add to the constraint 9 
extra-conditions MIN<t>, < <t>X, < [<t>,], where / = 1=9, then use the extreme condition for the equivalent error function 
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as 


in F = min ] - <I>X ()/ 0 - 5 ( < * > i + 4> x ,)] [< £ = 10 


to check the existence of valid solutions. In order to find the 


solution, it has to use the minimization condition of the equivalent penalty-function equal to zero or mnl ® 


together with 


F„ (T) = YPenahy(f(t)) + y) Penalty 

i=1 ;=1,5,6.9 


*;(*)-[*;] 

0-5 -|* y (t) + [*j]\ 


. From this, after filtration there are 14 Pareto 


optimal solutions, which are included in Table 7. In this scenario, the solutions are also represented visually by the chart in 


, <3> Sj -MIN<3> 

<3>/ i->- 1! — 

Figure 7a. The functions are converted into a relative form (unitless) ' ‘ , i = 1..9 ,j = 22..35, that they 

have the same scale in the chart. It shows that there is a similarity on magnitude among the functions 1, 2, 3, 5, 6, 7, 8, 9. 
Regarding the function 4, 3 solutions (S24), (S27) va (S33) are superior to the rest 11 solutions from the other 8 functions. 
Therefore, it points out that in order to meet conditions in scenario 3, it needs to opt for the solutions (S24), (S27) or (S33). 
The valid parameter domain in scenario 3, as shown in Figure 7b, is generally quite large, only for the third parameter the 
domain is rather narrow. 



-QCMAX QCMIN S.O ~ Sll-S12 SI3 SI4 S15-S16-S17-S18-S19-S20| 

- S2I __| 


Figure 6a: Intuitive Chart of Distribution of 12 Pareto Solutions 
with 9 Criteria of Saw Machine in Scenario 2 


0.0014517 _ 0.00147 


Figure 6b: Intuitive Chart of Solution Distribution in 
Accordance with Parameters - Scenario 2 


Table 7: Pareto Optimal Solutions in Scenario 3 


NO 

i 

1 

2 

3 

4 

5 

6 

7 

8 

9 

(S22) 

a-i 

0.03 

0.092196 

0.0014795 

0.06162 

0.17129 

0.29632 

1271.7 

2995.6 

- 

d>, 

1.0767 

0.99758 

-80.403 

-8608.4 

1271.7 

-2995.6 

-1737.4 

-55427 

0.0014795 

(S23) 

GC, 

0.034333 

0.065424 

0.0015012 

0.05336 

0.16591 

0.22663 

1296.5 

2848.4 

- 

ct, 

0.80846 

0.99549 

-80.647 

-7777.5 

1296.5 

-2848.4 

-1418.2 

-55089 

0.0015012 
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(S24) 

O-i 

0.033714 

0.076008 

0.0014981 

0.05627 

0.16574 

0.28538 

1278.1 

3000 

- 

ch, 

0.98148 

0.99391 

-80.144 

2.89E78 

1278.1 

-3000 

-1536.2 

-55027 

0.0014981 

(S25) 

a, 

0.034056 

0.068089 

0.0014996 

0.0544 

0.1573 

0.23564 

1290.5 

2988.1 

- 

ch, 

0.84012 

0.97771 

-80.382 

-5607.3 

1290.5 

-2988.1 

-1450.4 

-55278 

0.0014996 

(S26) 

a, 

0.033063 

0.068089 

0.0014998 

0.0544 

0.1573 

0.23564 

1290.5 

2988.1 

- 

ch, 

0.83858 

0.97573 

-81.205 

-5877 

1290.5 

-2988.1 

-1458.9 

-56064 

0.0014998 

(S27) 

a, 

0.032883 

0.079165 

0.0014975 

0.05813 

0.1507 

0.31068 

1267.9 

2861.7 

- 


1.0474 

0.96217 

-80.046 

2.89E78 

1267.9 

-2861.7 

-1584.1 

-55498 

0.0014975 

(S28) 

a, 

0.03 

0.089586 

0.001497 

0.06053 

0.14733 

0.29889 

1195.4 

3000 

- 


1.0737 

0.94966 

-82.014 

-5092.2 

1195.4 

-3000 

-1687 

-55024 

0.001497 

(S29) 

a, 

0.031363 

0.06151 

0.001494 

0.05499 

0.11795 

0.27034 

1294 

2870.6 

- 


0.86878 

0.89363 

-80.693 

-5036.8 

1294 

-2870.6 

-1409.6 

-57498 

0.001494 

(S30) 

a, 

0.031312 

0.074474 

0.0014887 

0.05695 

0.15672 

0.26807 

1292.4 

2998 

- 

ch, 

0.93188 

0.97106 

-80.943 

-31294 

1292.4 

-2998 

-1544.5 

-57129 

0.0014887 

(S31) 

a, 

0.030679 

0.076114 

0.0015002 

0.06036 

0.1403 

0.30555 

1145.2 

2822.8 

- 

ch, 

1.0176 

0.93696 

-81.298 

-9060.3 

1145.2 

-2822.8 

-1524.9 

-55329 

0.0015002 

(S32) 

a, 

0.030001 

0.078932 

0.0014938 

0.05974 

0.17046 

0.24249 

1193.3 

2907.8 

- 

ct>,- 

0.90344 

0.99592 

-81.744 

-5981 

1193.3 

-2907.8 

-1572.6 

-56244 

0.0014938 

(S33) 

a, 

0.032417 

0.079009 

0.0014832 

0.05606 

0.14229 

0.33038 

1298.2 

2806.4 

- 


1.0811 

0.94441 

-80.486 

2.89E78 

1298.2 

-2806.4 

-1562.1 

-55258 

0.0014832 

(S34) 

a, 

0.033288 

0.082257 

0.0014999 

0.05804 

0.16687 

0.28465 

1290.4 

2801.4 

- 

ch, 

1.0135 

0.99532 

-80 

-18022 

1290.4 

-2801.4 

-1625.7 

-55362 

0.0014999 

(S35) 

a, 

0.030766 

0.073734 

0.0014971 

0.06027 

0.17029 

0.30496 

1256.2 

2872.9 

- 


1.003 

0.99711 

-80.144 

-5127.6 

1256.2 

-2872.9 

-1565.8 

-58113 

0.0014971 



Figure 7a: Intuitive Chart of Distribution of 14 Pareto Solutions 
with 9 criteria of Saw Machine in Scenario 3 
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Figure 7b: Intuitive Chart of Solution Distribution in 



Accordance with Parameters - Scenario 3 
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Finally, there are 9 Pareto solutions available such as (S8), (S9), (S12), (S13), (S15), (S18), (S24), (S27), (S33) 
for all of three scenarios, as shown in Figure 8. Every solution corresponds to a particular manufacturing circumstance, 
defining the global behavior of the saw machine. Indeed, they are hardness, stability of saw-module, free vibration 
frequency, mass, geometry and technology. Based on these solutions, designer process used for the new type of frame saw 
machine can be carried out tluently among the stages such as detailed component design, verification of durability, 
hardness and stability, manufacture, installation, operation, ect. Yet, the designer can carry on opting for 3 representative 
solutions for three scenarios respectively in order to introduce them into manufacturing the saw machine. Those are the 
solutions (S9) - scenario 1, (S13) - scenario 2 and (S27) - scenario 3, design sketch of which is illustrated in Figure 9. 
Each parametric set presents a schematic diagram of the saw-module with the most important information such as 
sawblade dimension, part relative position, counterbalance position, etc. The designer can continue with the optimization 
problem of specific parts such as eccentric disks, housing parts and so on to produce the final drawing. However at these 
stages, the optimization design problem can be carried out easily by using the current available methods rather than VIAM, 
because it is mostly a single-objective optimization. 



Figure 8: Intuitive Chart of Distribution of 9 Pareto Solutions with 9 
Criteria of Saw Machine for all of Three Scenarios 



b) Solution (S13) c) Solution (S27) 


Figure 9: Illustration of Saw-Module corresponding to the Pareto Solutions 

CONCLUSIONS 

A multi-objective optimization problem including 8 control parameters, 9 functional constraints and 9 objective 
functions was elucidated properly by means of a visual interactive analysis method (VIAM) with an application of 
single-objective optimization techniques. Indeed, there were 28 Pareto optimal solutions, which have been determined for 
three manufacturing scenarios; hence the designer was able to make a suitable decision at every scenario. Also thank to 
VIAM, a new type of frame saw machine was designed successfully without inconsistencies among stages in design 
process such as concept, analysis, technology and multi-criteria decision-making processes. Besides, VIAM would 
definitely be reliable for designing another mechanical system apart from frame saw machine. 
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